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Abstract 

Dyson's integration theorem is widely used in the computation of eigenvalue correlation func- 
tions in Random Matrix Theory Here we focus on the variant of the theorem for determinants, 
relevant for the unitary ensembles with Dyson index (3 = 2. We derive a formula reducing the 
(n — fc)-fold integral of an n x n determinant of a kernel of two sets of arbitrary functions to a 
determinant of size k x k. Our generalisation allows for sets of functions that are not orthogonal or 
bi-orthogonal with respect to the integration measure. In the special case of orthogonal functions 
Dyson's theorem is recovered. 



1 Motivation 



Random Matrix Theory (RMT) has many applications in all areas of Physics and beyond (see e.g. 
the introduction of pQ). For the class of invariant RMT Dyson's integration theorem is at the heart 
of the method of orthogonal polynomials when computing all eigenvalue correlation functions exactly, 
for finite n x n matrices. The resulting expressions are then amenable to the large- n limit, in which 
universal RMT predictions follow. In the following we restrict ourselves to the integration theorem for 
determinants. Before presenting our generalisation thereof we briefly recall how it reveals all eigenvalue 
correlations in the unitary ensembles. 

We start by stating Dyson's integration theorem, as cited in pQ (Theorem 5.1.4). Given K(x,y) is 
a real valued function satisfying the following self-contraction property: 

dy K(x, y)K(y, z) = K(x, z) , 

^ dyK(y,y)=c. (1.1) 

Then it holds^l 

/dx\ det [K(xi,xj)] = (c — n+1) det [K(xi,Xj)] , (1-2) 
1 < i , j < n 2<i,j<n 

thus reducing the size of the determinant by one through the integration. The theorem also holds 
for kernels of orthogonal polynomials in the complex plane or for bi-orthogonal polynomials. A 
similar statement is true for quaternion valued kernels with the determinant replaced by a quaternion 
determinant (or Pfafnan). We refer to pQ for details as we will only consider the ordinary determinant 
case here. Iterating the integration theorem the following holds for an (n — fc)-fold integral: 

/n—k 

TT dxi det \K(xi,Xj)] = (c - n + 1) . . . (c - k) det [K(xi,Xj)} . (1.3) 

- LJ - Ki,j<n n—k+Ki,j<n 

1=1 ~ ~ J - 

It is this form that we will generalise as it is most useful when computing correlation functions in 
RMT. We emphasise that on the right hand side (rhs) the determinant has reduced to size k x k over 
the same kernel. 

The application of eq. (|1.3j) to the unitary ensembles goes as follows. Suppose we have a set of 
orthogonal polynomials Pk{x) = x k + 0(x k ~ 1 ) of order k in monic normalisation satisfying 



/ 



dx w(x)p k (x)pj{x) = h k 5 kj . (1.4) 



Here w(x) is a positive weight function such that all moments exist. From the polynomials we can 
construct orthonormal wave functions 

<p k (x) = ht 1/2 w(x)V 2 p k (x) (1.5) 

and the following kernel 

n-l 

K n (x,y) = ^2ipj(x)ipj(y). (1.6) 

j=Q 



1 The symmetry property K{x,y) = K(y,x) stated in [T] is not necessary, as can be seen from the proof there. 



2 



It satisfies Dyson's theorem above with c = n. In the symmetry class of unitary invariant RMT, the 
partition function is given in terms of the joint probability distribution (jpdf ) of eigenvalues as 

/n 

~[dxi w(xi) A n (x) 2 . (1.7) 
i=l 

The Vandermonde determinant in the integrand, 

A n (x)= det [aJ- 1 ] = det bi-l(^)] , (1-8) 



l<i ,j<n 



l<i,j<n 



can be replaced by a determinant over an arbitrary set of monic polynomials. If we choose the 
orthogonal ones we can rewrite the jpdf and thus the partition function after few manipulations as 

/n 
Idxi hi-i det [K n (xj,x k )}. (1.9) 
■ LJ - i<j,k<n 
i=i ~ J ~ 

It immediately follows from Dyson's theorem that Z n = n\ Y\7 = x hi-i- Moreover, following eq. (|1.9p all 
/c-point eigenvalue correlation functions given by n — k integrations over the jpdf can be immediately 
read off: 

1 f n 

Rk{xi, . . . ,Xk) = 7 ttt / dxi det [K n (xi,Xj)] = det [K n (xi,Xj)] . (1.10) 

n — K ! / - LJ - Ky<n l<i,j<k 

v 1 J i=k+l 

In the large-n limit the size of the determinant on the rhs remains fixed, and the kernel can be easily 
evaluated using the Christofiel-Darboux identity for orthogonal polynomials on R. 

Since in this example the choice of orthogonal polynomials was entirely at our disposal, why should 
we choose polynomials that are not orthogonal with respect to the weight function, or the integration 
range? The reason is that we are not always able to choose the polynomials to be orthogonal. One 
example where such a situation occurs is in the Schwinger model [2]. A second example, being in 
a different symmetry class, appears when considering the Ginibre ensemble with real non-symmetric 
matrices [3] . Integrating out all real eigenvalues one arrives at the Pfaffian of the so-called D-kernel of 
the Gaussian Orthogonal Ensemble [I], integrated over a non-Gaussian weight function in the complex 
plane. Consequently the self-contracting property eq. (jl.ip is not satisfied then. 

For this reason we propose a generalisation of Dyson's theorem for determinants without imposing 
any orthogonality condition, and we restrict ourselves to real integrals for simplicity. The generalisation 
to integrals over C is straightforward. A counterpart for integrating Pfaffians of a non self-contracting 
bilinear has been proved in [3] in the special case when all variables are integrated out. 

2 Results 

Let each {4>j{p)} and {ipj(q)}, j = 1, . . . , n be a set of linearly independent, integrable functions, such 
that all integrals J dx <f)i(x)ipj(x) exist. For these two sets we define the following bilinear function: 



Qn(x,y)= ^2</>j(x)il>j(y) . (2.1) 



Then the following holds: 



j(x) may or may not be linear combinations of the functions tpj(x) 
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Theorem 1: 



n I ]I dx l ^det [QniPhQj)} 

O / Ki,j<n 

; =1 



(n-/c)! det <?j)] , & = 0, 

n— fc+l<i j'<n 



(2.2) 



where we have set pi = qi = Xi, i = 1, . . . ,n — k for all integration variables. 
The kernel IC n (p, q) on the rhs is given by 



1 n 



a=l 



J dx 4>i(x)i/ji(x) 
f dx 4> n {x)ipi(x) 



4>l (p)V'a (?) ■■■ J dx 4>i (x)lp n (x) 
<Pn(j>)lpa(q) ■■■ f dx <j) n (x)lp n (x) 



(2.3) 



where the sum runs over the a-th column replacing the integrated functions by unintegrated ones. 
The normalisation C on the left hand side (lhs) is given by 



C = det 



J dx cpi(x)^i(x) • 
/ dx 4> n (x)ipi(x) 



■ f dx (p 1 (x)ip n (x) 
J dx cj) n (x)ijj n (x) 



(2.4) 



Thus we have reduced an (n — £;)-fold integral over an n x n determinant to a k x k determinant of a 
single kernel, consisting of a sum of n x n determinants containing only single integrals. 

The rhs of our Theorem 1 can be interpreted as a generalised kernel having 2k variables. If we 
define 

Xpu ■■■ ,Pk',Qh ■■■,%) = 77 det [Knipuqj)] , (2.5) 

K\ l<i,j<k 

these satisfy the following generalised self-contraction property (see eq. (jl.ip ): 
Theorem 2: 

J dqi...dq k K^\pi,...,p k \qi,...,q k )}C^\qi,...,q k -ri,...,r k ) = K$\pu 

h---dq k 1C ( £ ) (qi,...,q k -,qi,...,q k ) = 



,Pk',n, 



,r k ) 



!■ 



(2.6) 



In particular the kernel defined in eq. (|2.3p . /ci fc_1 ^ (p; q) = K n {p;q), is self-contracting. 

Let us make a few remarks. First, the bilinear Q n (x,y) of the set of functions is in general 
different from the kernel on the rhs: Q(x,y) ^ IC n (x,y). In particular it is not self contractive in 
general: J dy Q(x,y)Q(y, z) ^ Q(x,z). In the case of orthogonal functions, J dx (j) k (x)ijjj(x) = 5 k j, 
we obviously get back Q n (x,y) = IC n (x,y). Then Dyson's theorem applies, as in the example in the 
previous section. 

Special cases of Theorem 1 were previously known. For k = it goes back to C. Andreief in 1883 
as cited in [1], after multiplying with the normalisation C: 



T dx\ det [4>j(xi)] det [ipj(xi)] = n\ det 

Ki,j<n l<i,j<n Ki,j<n 

1=1 



dx 4>i(x)tpj(x) 



(2.7) 



The identity for k = 1 was stated and used in [2] but no explicit proof was given. Furthermore let 
us point out that for k = n there are no integrations, thus equating the determinant of the bilinear 
function and of the kernel. 

All statements above also hold when the normalisation eq. (|2.4p accidentally vanishes, C = 0, as 
will be indicated below. This cannot happen for Dyson's integration theorem. 
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3 Proofs 



The proof of Theorem 1 will go in three steps, taking the known result for k = for granted. In step 
i) we prove the Theorem for k = 1, relating to the definition (12, 3D . In step ii) we show that the kernel 
JC n (p,q) satisfies the self-contraction property eq. (|2.6p . Theorem 2 for k = 1. In the last step iii) we 
prove Theorem 1 for k = n without integrations. Because of the self-contraction property of lC n (p,q) 
we can then apply Dyson's theorem to the rhs to show all the remaining cases. Theorem 2 for k > 2 
will then be shown in the second part. 



as follows: 



step i): To prove k = 1, in a first trivial step we can replace the determinant of the bilinear function 

n 

Y<t>a{Pi)ll>a{qj 



det [Q n (Pi,qj)] 



det 

l<i ,j<n 



,o=l 



det [0 a (ft)] det [V>a(<7j)] • 

l<a,j<n lS a j£ n 



(3.8) 



Inserting this into the lhs we can expand both determinants with respect to the last, unintegrated 
column: 



.. ,. n— 1 -. n— 1 / n \ / n ^ 

C Hdxi drt [^(pi)] det [^fe)] = ~ / 11^™ E^CP") ? E^(</nX? 

J=l - a > 4 - n _a,J_n J m=1 I , =1 I \l=l / 



(n-1)! 



C 



1,3=1 

= (n- 1)1 JC n (p n ,q n ) , 
where p% = q% = Xj for i = 1, . . . , n — 1. We have introduced the minors 

C/ = (-) n - 1+J det[^(x fc )] and cf = (-) n " 1+i de^a*)] 



(3.9) 
(3.10) 



These contain n — 1 functions each, and all variables x±, . . . , x n -\ are integrated. Thus for each product 



Cfcf we can apply the formula for A; = by C. Andreief eq. (|2.7p . with the resulting minor 



Cu 



det 



dx 4>i(x)ip k {x) 



(3.11) 



In the last step the sum in eq. (|3.9p can be precisely written as the sum over determinants in eq. 
(12. 3p . each expanded with respect to the a-th column. 



step ii): To derive the self-contraction property for JC n (p,q) we simply insert the definition eq. (|2.3p . 

<Pi(p)ipa{q) ■■ 



applying the short hand notation < i,k >= J dx 4>i(x)^k(x): 
q JC n {p,q)K. n (q,r) 



1 

C2 



n 


< 1, 


1 > 


E det 






0=1 


< n 


1 > 


n 


' < 1, 


1 > 


E det 






6=1 


< n 


1 > 




TT2 E det 



< 1,1 > • 

< n, 1 > ■ 



4>i{q)A{r) 
4>n{q)A{r) 




< n, n > 
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Here we simply observe that in each product of determinants the common factors ipa(q) and ipb(^) 
can be taken out of the columns a and b respectively, and can then be multiplied into the columns of 
the other determinant. In a second step the integral J dq can now be taken inside the 6-th column of 
the second determinant containing <pj{q)ip a {q), resulting into < j,a >. This leads to a column already 
present and thus a vanishing determinant, unless we have a = b. The resulting normalisation cancels 
one power of C to reproduce the kernel. 

step iii): We prove Theorem 1 for k = n when all integrations are absent. It is easily seen when 



expanding the kernel inside the determinant on the rhs, using the formula in the last line of eq. (|3.9|) : 

1 n 



det [K,{pi,qj)\ 

±<t,j <n 



det 

l<i j<n 



Lk=l 



det [(f> k (pi)] det [ipi(pj)] det — Q fc 

l<t,k<n 1<'J<"- l<(,fc<n Cv 



77 , det [Qn(Pi,qj)} ■ 



(3.13) 



In the first step we used that the determinant of the matrix product is the product of the determinants. 
Furthermore, the minors are just the matrix elements of the inverse matrix, Cik/C = (C _1 );fc, and we 
used eq. (13. 8j) . Because of ii) we can now apply Dyson's integration theorem to the determinant of 
the self-contracting kernel fC n {pi, qj). Using its normalisation, J dq IC n (q, q) = n which can be trivially 
seen, we arrive at Theorem 1 for all h = 0, . . . , n. 

Finally we deal with the case C = 0. The Andreief formula eq. f)2.T[) still holds and vanishes 
identically. To make Theorem 1 nonsingular we multiply it with C k . Step i) proving k = 1 remains 
true and in general nonzero, using eq. (|2.7|) for size n — 1. For k > 2 the lhs of Theorem 1 times C k is 
zero. If the matrix under the determinant C eq. (|2.4p has rank n — 2 or less, all determinants inside 
the kernel eq. (12. 3D vanish and make the rhs identically zero. For rank n — 1 it is easy to see that the 
kernel eq. (|2.3|) factorises into two functions of p and q, thus having a vanishing determinant on the rhs. 



It remains to prove Theorem 2 for k > 2. The proof goes as follows. On the rhs of f|2.6|) we 

rk\ 

substitute each of the two kernels K/ n as a determinant of single kernels from their definition (12. 5|) . 

Using the standard representation of a determinant, we get: 

k 

fcn(Pa(i),qi)fcn(qi,r a ,(i)) 



ihs= w/^-^^ ( - i)CT+,T 'n 



(3.14) 



The fact that IC n (x,y) is self-contractive allows us to do all the integrals over to obtain 

k 

lhs = W^^ 1 ^' IlM*M0>Mi)) = h,M c J C ^ T j) = K ( nHpu-,Pk;ri,...,r k ). (3.15) 



(k\) 



(k) 

The normalisation of the generalised kernel ICn follows directly from the normalisation of the single 
kernel )C n (p,q) which is n, and the repeated application of Dyson's integration theorem to a k x k 
determinant integrated k times. This way we get the binomial coefficient: 



J dqi...dq k ld£\qi, . . . , q k ; qi, 



,qk) 



(n-k)\ k\ 

The same arguments as at the end of the previous proof apply for C = 0. 



(3.16) 
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4 Conclusions 



We have shown how to reduce any number of integrations over a determinant of a bilinear function of 
non-orthogonal functions to a smaller determinant of a self-contracting kernel containing only single 
integrals. This makes the large-n limit feasible in such a general setting at least in principle, given 
the single integrals can be evaluated. Due to the fermionic nature of the Vandermonde determinant 
other applications than the mentioned Schwinger model should exist. Our result gives hope that an 
analogous Pfafnan integration theorem with some variables unintegrated also exists. 

After writing up this paper we learned from P. Forrester that the first part of our result was derived 
independently by Rains [5] in the context of symmetrised increasing subsequences. His alternative 
proof is formulated in terms of the Pfafnan of an antisymmetric matrix kernel. In contrast, our proof 
illuminates the close relation to Dyson's theorem. 
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